An exact mapping between the states of arbitrary N-level quantum systems 
and the positions of classical coupled oscillators 
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The dynamics of states representing arbitrary N-level quantum systems, including dissipative 
systems, can be modeled exactly by the dynamics of classical coupled oscillators. There is a direct 
one-to-one correspondence between the quantum states and the positions of the oscillators. Quantum 
coherence, expectation values, and measurement probabilities for system observables can therefore 
be realized from the corresponding classical states. The time evolution of an N-level system is 
represented as the rotation of a real state vector in hyperspace, as previously known for density 
matrix states but generalized here to Schrodinger states. A single rotor in n-dimensions is then 
mapped directly to n oscillators in one physical dimension. The number of oscillators needed to 
represent N-level systems scales linearly with N for Schrodinger states, in contrast to N 2 for the 
density matrix formalism. Although the well-known equivalence (SU(2), SO(3) homomorphism) of 
2-level quantum dynamics to a rotation in real, physical space cannot be generalized to arbitrary 
N-level systems, representing quantum dynamics by a system of coupled harmonic oscillators in 
one physical dimension is general for any N. Values for the classical coupling constants are readily 
obtained from the system Hamiltonian, allowing construction of classical mechanical systems that 
can provide visual insight into the dynamics of abstract quantum systems as well as a metric for 
characterizing the interface between quantum and classical mechanics. 



I. INTRODUCTION. 

The density matrix formalism [J-Q provides a straight- 
forward procedure for predicting quantum dynamics. At 
any given time, the density matrix provides a complete 
statistical characterization of the system in terms of the 
mean values of measurable states. It includes both the 
quantum uncertainty in predicting the results of single 
measurements on pure states (even though such states 
represent maximal possible information for the system) 
and the classical uncertainty in measurements on mixed 
states of less than maximal information. 

Although the theory needs no supporting visual model 
for its application, the Liouville-von Neumann equation 
governing the time evolution of the density matrix pro- 
vides little physical insight into system dynamics. There 
has therefore been considerable effort towards represent- 
ing, where possible, quantum systems using more intu- 
itive classical models. Of particular influence and im- 
portance is the classical representation for quantum two- 
level systems [H , sometimes referred to as the Feynman- 
Vernon-Hellwarth (FVH) Theorem. The behavior of any 
quantum mechanical two-level system can be modeled by 
classical torque equations, providing a one-to-one corre- 
spondence between the time evolution of the system and 
the dynamics of a spinning top in a constant gravitational 
field or a magnetic moment in a magnetic field. 

Work by Fano 3] published concurrently with the FVH 
result also provides geometrical interpretation of spin dy- 
namics for more complex quantum systems. The density 
matrix for an N-level system is represented as an expan- 
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sion in Hcrmitian operators, resulting in a vector with 
real components. The time development of this gener- 
alized Bloch vector is a real rotation in a hyperspace of 
(N 2 — 1) dimensions. Constants of the motion can be de- 
rived [Bllll that constrain the system's dynamics and pro- 
vide physical insight. However, the states of the system 
as given by the components of this vector (also referred 
to as a coherence vector Q and, more recently, a Stoke's 
tensor Q), do not evolve in a physical space amenable to 
visualization, with its attendant advantages, except for 
the case N — 2. 

Thus, no completely general mapping has been realized 
that yields a one-to-one correspondence, similar to the 
FVH result, between the states of a quantum-mechanical 
N-level system and classical dynamical variables, provid- 
ing the possibility for direct mechanical insight into the 
dynamics of abstract quantum systems. Analogies be- 
tween quantum and classical systems have been noted 
[§-r20j almost from the beginning. But exact equivalence 
between the quantum and classical equations of motion 
has been obtained only for certain limiting conditions 
[2l| - |28| such as weak perturbations of the system (weak 
coupling limit) and the aforementioned 2-level systems. 

Recently, the possibility of exact representations of N- 
level quantum systems in terms of classical coupled oscil- 
lators, with no restriction to weak coupling, was demon- 
strated [2{|. However, the formalism is limited to real, 
invertible Hamiltonians applied only to pure states. It 
is insufficiently general to represent statistical mixtures 
and density matrix evolution classically, as well as open, 
dissipative systems. 

In the present work, a very simple approach is pre- 
sented for mapping an arbitrary N-level quantum system 
to a system of coupled harmonic oscillators, available for 
over a decade. It both complements and augments the 
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results presented in [29j | and is made more relevant by 
that study. The salient features representing the dynam- 
ics of N-level systems as rotations in Liouville space are 
reviewed first. The desired one-to-one correspondence 
between the states of the quantum system, represented 
as a density matrix, and classical dynamical variables 
is provided by a mapping representing harmonic oscilla- 
tors. The quantum states, either pure or mixed, are rep- 
resented exactly by the time-dependent displacements of 
classical coupled oscillators. 

Hilbert space rotations are reviewed next to general- 
ize this option for representing quantum spin dynamics 
classically. Results inspired by [29| are presented with 
no restriction to real, invertible Hamiltonians. An exact 
mapping of Schrodinger states to the physical displace- 
ments of coupled oscillators is provided, in contrast to 
the displacement and velocity originally required. This 
approach to representing spin dynamics is then extended 
to mixed states. Whereas N 2 — 1 oscillators are needed 
most generally to represent density matrix dynamics for 
an N-level system classically, Schrodinger states require 
at most 2N oscillators. 

Open (dissipative) N-level systems are considered next, 
showing they also can be exactly and simply modeled as 
classical coupled oscillators. The present treatment re- 
veals the necessity for negative couplings in closed sys- 
tems as well as antisymmetric couplings in open systems. 

The paper closes with illustrative examples of the 
quantum-classical mapping. 



II. TIME EVOLUTION OF QUANTUM 
N-LEVEL SYSTEMS 

A brief synopsis of formalisms for representing the dy- 
namics of N-level systems is presented. The standard Li- 
ouville space and Hilbert space representations are con- 
sidered first, enabling simple generalizations that lend 
themselves to a classical interpretation. In all cases, the 
time evolution of the system can be reduced to the form 

S(i) = tf(t)*(0). (1) 

The representation chosen determines the particular 
forms for $ and the propagator U{t). For notational con- 
venience and interchangeability of energy and frequency 
units, h is set equal to 1 in what follows. 



which defines U(t) — e" lHt . 

The time evolution can be related to a rotation by first 
expanding p in terms of a complete set of basis operators 
3] . Orthogonal bases are particularly convenient and are 
typically normalized for further convenience. Denoting 
the basis elements as for state i and requiring only 
that the basis be orthonormal gives 



{e i \e J )= Tr(e| e,-) = 



(4) 



where the inner product for the vector space comprised 
of matrices is given by the operator Tr, which returns the 
trace (sum of diagonal elements) of its argument, and f 
denotes the operation of Hcrmitian conjugation. In lieu 
of explicitly normalizing the e^, the inner product can 
be defined with the appropriate factor multiplying Tr. 
Summing over repeated indices gives p as 

P = r j£j (5) 
where the coefficients in the expansion are the projection 



onto the basis states. Each 



Eq. ©, 



r j = (e j \p)=Tr(etp), 

is thus the expectation value of the quantum state ej . 
Then 



(6) 



Tr(eJ p) 



!Tr(ej[tf,p]) 



(7) 



Expanding the commutator, using Tr(AB) = Tv(BA) = 

2 J J ! - ; 



[Tr(AB)t]* and [ ej , e 3 ]t = [e] , ] gives 



= -iTr(el[H,e 3 ]) 
= iTr([ele j ]H) 
= -{ l Tr([e],e l }H)y 



(8) 



in terms of its complex conjugate elements, denoted by *. 
Thus, fl = — fJt is antihcrmitian and can be diagonalized. 
The evolution of the density matrix is given by 



r = fir, 



(9) 



with solution 



A. Liouville equation 

The Liouville- von Neumann equation for the time evo- 
lution of a density matrix p governed by system Hamil- 
tonian H is 

p=-i[H,p], (2) 

with formal solution 

p{t) = e- lHt p(0)e zHt 

= Up{Q)U\ (3) 



r{t) = e ut r(0) 



(10) 



The propagator U(t) = e~ ot , and therefore W = U^ 1 
is unitary, since & = —ft. Thus, Eq. (j9]) represents a 
rotation, albeit still most generally in complex space. 



1. Rotation in real space 

In the case where the orthonormal basis states are Her- 
mitian operators, the components of the density matrix 
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(coherence vector or Stokes tensor, in this case) are real, 
and f2 is also a real antisymmetric matrix. The rotation 
of Eq. ([9]) is then a rotation in real, multidimensional 
Cartesian space, which is the generalization [1, H| to N- 
level systems of the result obtained in Q]. The quan- 
tum dynamics are thus fully classical in the additional 
dimensions exceeding 3D physical space. However, clas- 
sical rotations in more than three dimensions are only 
marginally less abstract than rotations in Hilbert space 
and Liouville space. More accessible insight can be ob- 
tained by mapping the real-valued quantum states to 
physical space as follows. 



2. Exact mapping to classical coupled oscillators 

A textbook exercise for deriving the Larmor precession 
of a spin-1/2 in a static magnetic field B differentiates 
the first-order derivative in Ehrenfest's theorem to obtain 
a harmonic oscillator equation for the time evolution of 
expectation values of the spin components transverse to 
Bq. Similarly, differentiate Eq. (|9|) again to obtain 



r = il 2 r. 



(11) 



Since J7 2 is real, symmetric, and diagonalizable, the so- 
lution is readily written in terms of the normal-mode so- 
lutions obtained by diagonalizing Q 2 . The eigenvalues of 
antihermitian Q are pure imaginary, resulting in negative 
eigenvalues — lu 2 (a = 1, 2, 3, . . . , n) for n x n matrix Q, 2 . 
The n distinct eigenvectors | w a ) constitute a basis set 
satisfying the completeness relation ^ a \ uj a ) ( uj a \ = 1 
(the identity element). In this eigenbasis, Eq. (fTTj) for 
each component r a — {uj a \r) becomes simply 



with standard harmonic oscillator solution 



r a (t) = r a (0)cosuj a t + 



r«(0) 



sin u a t 



(12) 



(13) 



and r dependent on r according to Eq. (JSJ, giving 



"(*)> = E Iw«>(wa|r(t)) 



= E |Wa)(Wa 
o=l 

= U(t)\r(0)). 



cos uj a t + Q 



Sill U! a t 



"(0)> 
(14) 



Representing the propagator in the eigenbasis where 
il is diagonal gives Eq. ([T2j) . Using the representation 
for non-diagonal fl and its eigenvectors to calculate U(t) 
gives the physical displacements r.i(t) for each of the n 
oscillators. This solution for U (t) must be identical to the 
propagator given in Eq. (|10j). It is included primarily for 
consistency in the presentation, but also to emphasize the 
fundamental differential equation under consideration is 



first order and only requires specification of r(0) (or r 
evaluated at any other fixed time). 

To complete the explicit identification of Eq. (jlll) with 
mechanical oscillators, consider equal masses, m, on a 
frictionless surface, with mass connected by spring of 
stiffness kij = kji to mass rrij — 1, 2, 3, . . . , n), as in 
Fig. Q] for an illustrative case n = 3. The classical matrix 

ki3 



*12 *23 



r 2 



FIG. 1: Schematic of three masses at equilibrium positions 
n = coupled with springs of stiffness kij. 

f2ci relating the displacement from equilibrium of the i th 
mass to its acceleration, as in Eq. (fTTj) . is 



i=i 



i = i 



(15) 



The above expression assumes reciprocal couplings 
kij = kji and employs the absolute value in the sum 
to accommodate negative couplings. A system of pendu- 
lums consisting of masses attached to rigid rods can be 
coupled negatively by attaching a spring to rod i below 
the fulcrum of oscillation and to rod j above the fulcrum. 
Displacing mass rrii to the right exerts a force on rrij to 
the left, ie, the coupling kij < 0. A pendulum can be 
inverted with its mass above the fulcrum to implement 
ku < 0. Inverting transformers can be used to implement 
negative couplings in LC circuits. 

Setting Slci = fl 2 gives the spring constants 



m 



i 7^ 3 
i =3 



(16) 



in terms of the matrix 51 (squared) representing the quan- 
tum system, as derived from Eq. There is thus a one- 
to-one mapping of the quantum states to the oscillator 
displacements embodied in both systems in 7"j(t). Given 
the initial states r.j(0) of the system, the necessary fi(0) 
follow from Eq. ([9]). 

This mapping is very general. It is not limited to 
particular values of the spin, numbers of interacting 
spins, specific forms of the commutation relations, or rel- 
ative fractions of mixed and pure states comprising p. 
An N x N density matrix generates N 2 components in 
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Eq. ([5]), which requires N 2 oscillators. The static com- 
ponent of the identity element can be eliminated, and 
the structure of the Hamiltonian may generate evolution 
only in a smaller subspace of states, further reducing the 
number of required oscillators. 

For pure states, the time-dependent elements Cj(i) 
comprising the state vector can readily be obtained, if 
desired, from p reconstructed in matrix form using the 
Ti{t) and Eq. ([5]). Each resulting element pij is equal to 
CiC*. Assigning any one of the q to the square root of 
Pa sets the arbitrary global phase of the pure-state ele- 
ments. In terms of this real Cj, the remaining Cj are equal 

to pji/Ci- 

B. Schrodinger equation 

The solution 

|*(i)) = e- ifft |*(0)). (17) 
to the time-dependent Schrodinger equation 

im)) = Hm)) as) 

represents a rotation of ^(O) ) in Hilbert space, since H 
is Hermitian and the propagator U(t) — e~ tHt is uni- 
tary. Most generally, H and the components Ci of | VP ) 
in a chosen basis are complex, so there is no classical 
interpretation for the time evolution of the state. 

In [111, the authors derive a methodology for repre- 
senting the complex Cj in terms of the displacement and 
velocity of classical coupled oscillators for the special case 
of real, invertible %. These restrictions can be removed, 
as will be shown in what follows. 



1. Rotation in real space 

The results of the preceding sections can be extended 
to the Schrodinger equation using only the quantum 
Hamiltonian, in contrast to the approaches in y| HI, [29| 
which start with the real, classical Hamiltonian % and 
relate it to the quantum Hamiltonian H. The expression 
for TL in terms of H then represents classical coupled 
harmonic oscillators only for real H 29] . 

Instead, start with Eq. (fT8|) and represent the compo- 
nents of I ) as vector c, giving 

c(t) = -iHc(t). (19) 

Write complex c = q + ip and complex H = Q + iP. For 
real H, q and p are conjugate variables [29( but are not 
most generally so. Equating the real and imaginary parts 
after performing the multiplications in Eq. (|19[) recasts 
the Schrodinger equation as a real equation 

(l) = (-Q p) (?) 

= n(j) (20) 



in the form of Eq. ([9]) with Qt — _ 57 again antiherme- 
tian, since Hermitian H requires = Q and P^ = —P. 
The results and implications for real rotations then follow 
from § 111 A II and the discussion following Eq. ^ . 

2. Exact mapping to classical coupled oscillators 
Similarly, differentiating Eq. (|2T))) gives 




in the form of Eq. (jlip for real symmetric (Hermitian) 
matrix f2 2 constructed from the real and imaginary parts 
of H 2 . The mapping of q and p to mechanical oscillators 
then follows from ^11 A 21 

For complex N x N Hamiltonian, there are thus most 
generally 2N mutually coupled oscillators. There can be 
fewer oscillators and no mutual coupling between spe- 
cific oscillators, depending on the structure of H. The 
displacements q and p provide an exact one-to-one map- 
ping to the real and imaginary components, respectively, 
of the quantum state ). The state c(0) = q(0) +ip(0) 
uniquely determines the initial displacements, with the 
initial velocities then given by Eq. (l20l) . 

The present treatment produces the results in [29l | in 
the case of real H = Q, P — 0. Note that q and p 
then evolve independently according to the same prop- 
agator, with no mechanical coupling between q and p 
oscillators. The initial conditions are the only difference 
in the solutions. Calculating p = PL~ x q separately as in 
[2^ | imposes an additional unnecessary restriction that 
H, already constrained in [2{| to be real, must be invert- 
ible (ie, no eigenvalues equal to zero). 

3. Extension to mixed states 

A statistical mixture can not be represented in terms of 
a state | ^ ) , but is written in terms of the probability pk 
for being in each of the possible states | ) as a density 
matrix 

P(t) = $>*!**(*)><**(*) I (22) 
k 

which evolves according to Eq. (|2|). The results in [29| 
and extensions in the previous section are limited to pure 
states evolving according to the Schrodinger equation. 
The methodology would appear to be inapplicable to 
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mixed states. However, the density matrix representing 
a given system is not unique. It is an average over the 
TV constituents comprising a macroscopic system, which 
can be astronomically large, precluding an exact deter- 
mination of the exact state of each of the A/" constituents. 

But the identical density matrix can also be con- 
structed from a completely specified set of N -C Af non- 
interacting pure states, with the N 2 elements of p deter- 
mined from measurable macroscopic (average) properties 
of the system, such as energy or polarization. In that 
case, both the weights pk and corresponding states are 
known exactly, so each | \&fc ) can be used independently 
to construct a set of coupled oscillators representing the 
components cf \t) of | ^fe(t) ). In lieu of density matrix 
evolution pit") = Up(0)W , the simpler and more efficient 
Schrodinger evolution | *(<) ) = U\ ^(O) ) can be applied 
to each pure state | ) comprising p in Eq. (|22|) . Sub- 
sequently, the weights pk can be used to calculate expec- 
tation values, measurement probabilities, or reconstruct 
the density matrix at later times t if desired. 

In addition, as shown in [2fj, at least one of the | VPfc ) 
comprising the initial density matrix is redundant and 
can be removed from the calculation, since it provides 
a relatively uninteresting constant contribution to the 
system dynamics. Choose one of the weights, for ex- 
ample, pi. The density matrix can be rewritten as p\ 
times the identity element plus a "pseudo" density matrix 
constructed from the | "Jfc ) with weights (pu — Pi)- The 
term that is proportional to the identity element doesn't 
evolve in time under unitary transformations and can be 
ignored. 

Thus, the state | *i ) has been removed from the den- 
sity matrix, along with any other | vf^ ) that had original 
weights pk — Pi- In the general case of m > 1 degener- 
ate statistical weights pk, only m < N of the VPfc ) are 
required. The number of oscillators is correspondingly re- 
duced. Choosing the weight with the largest degeneracy 
provides the maximum reduction. The explicit contribu- 
tion of each pure state to the system dynamics is readily 
apparent in utilizing this approach. The density matrix 
at any given time is easily reconstructed as described in 



III. DISSIPATIVE SYSTEMS 

The modifications necessary to model open systems 
as a set of damped oscillators can be found very gener- 
ally, with minimum detail concerning the relaxation for- 
malism. The Wangsness-Bloch equation expressing the 
evolution of the density operator in the presence of re- 
laxation adds a relaxation operator term to Eq. ([2]) that 
operates on the density matrix 0, |3l| . Expanding p in 
a basis of orthonormal operators as in Eq. ([5]) gives the 
real equation 



which can be transformed to a homogeneous equation as 
in |3J . The relaxation matrix R must be symmetric for re- 
laxation elements that act symmetrically between states 
of the system, with diagonal elements providing auto- 
relaxation rates and off-diagonal elements giving cross- 
relaxation. The function F adds a constant term incorpo- 
rating the asymptotic decay of the system to the steady 
state in terms of the equilibrium state r eq . Without this 
term, the solution decays to zero. 
Differentiating again gives 

r = (n + R)r 

= + R)r + F}+ Rf, (24) 

ie, a set of coupled oscillators with a velocity-dependent 
friction term and a constant applied force flF. A con- 
stant force in the harmonic oscillator equation merely 
shifts the origin of the coordinates. However, the ma- 
trix multiplying r, which determines the mechanical cou- 
plings as given in Eq. (|16p , is no longer symmetric due to 
the sum of antisymmetric f2 and symmetric R, resulting 
in non-reciprocal off-diagonal couplings. 

The precise role of non-reciprocal couplings in a clas- 
sical model for quantum dissipative systems can be clar- 
ified by eliminating r to obtain 

r = (fi + R) 2 r + + R)F 

= r 2 r + TF, (25) 

a set of ideal (frictionless) coupled oscillators subjected 
to a constant applied force. In this case, however, the 
matrix T 2 is the sum of symmetric VI 2 + R 2 and anti- 
symmetric flR + Rtt. The former term corresponds to 
a set of undamped oscillators with symmetric couplings 
kij = kji f ^II A 2|) . modified in comparison to no relax- 
ation by inclusion of R 2 . The normal- mode frequencies 
are also modified accordingly. 

Damping is provided by the antisymmetric part of T 2 , 
which gives antisymmetric couplings 7^ = —"fji and total 
coupling Kij = + jij. The Yy therefore represent 
couplings connected in parallel with the symmetric kij 
and can be implemented, in principle, using magnetic 
materials and magnetic fields. For a given positive 7^, 
a positive displacement of mass m,j results in a positive 
force on m,-. The resulting positive displacement of rrii 
provides a negative force on rrij due to jji < which 
opposes the original displacement of rrij and damps the 
motion. Stated differently, energy transferred from rrij 
to rrii is not reciprocally transferred back from rrii to rrij , 
and the motion is quenched. 

The inhomogeneous term TF in Eq. (|2"5"|) can, equiv- 
alently, be included in an augmented matrix T formed 
by appending a column TF to the right of T 2 and then 
adding a correspondingly expanded row of zeros at the 
bottom. The vector r is then augmented by including a 
last element equal to one to obtain the equivalent homo- 
geneous equation 



r = + Rr + F(r eq ), 



(23) 



u _ ~ _ 
— r — r v 

dt 2 



(26) 
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The asymmetry of T generates unphysical couplings that 
are not a problem theoretically but preclude a real, 
physical model. However, T is readily written as the 
sum of symmetric T$ = (T + T^)/2 and antisymmetric 
Ta = (f — T')/2, which determine the symmetric cou- 
plings kij and antisymmetric couplings 7y, respectively, 
as above. 

In comparison, the Schrodinger equation can only in- 
clude relaxation in certain special cases amenable to com- 
plex energies in the Hamiltonian. A typical application 
is the coupling between stable and unstable states and 
the resulting lifetimes of the states. An example relat- 
ing velocity-dependent damping of classical oscillators to 
a Schrodinger equation treatment was provided in [29j 
in the weak-coupling limit. However, neither this ap- 
proximation nor the required complex energies can be 
applied more generally. Even a simple two-level system 
with relaxation dynamics described by the Bloch equa- 
tion cannot be addressed by the Schrodinger equation 
and requires the density matrix approach. 



U(t) = e -i"ot e -ivcr 

_ e -iu t [ cosu t _ iQi ■ cr smut] 

= e "" 0t ( -b* a* ) • ( 28 ) 

The parameters a, b obtained from expanding d> • cr and 
using the matrix forms for the cr, are 

a = cos cot — i c2?3 sin cot 

b = — {t02 + i&i) sinwi, (29) 

recognizable from classical mechanics as the Cayley-Klein 
parameters for a rotation by angle 2u>t about uj. 

Evolution of Schrodinger state | ^ ) <H> (ci, c-i) proceeds 
according to Eq. (|17[) . with the corresponding density 
matrix states pij — C{C* evolving according to Eq. ([3]). 
The equivalent classical evolution is considered next. 

B. Classical representation (Liouville equation) 



IV. ILLUSTRATIVE EXAMPLES 

Simple two-level systems are used as a prototype for 
implementing the quantum-classical mapping. Although 
they are already known to be representable by classical 
rotations in three-dimensional physical space, they pro- 
vide sufficient detail to clarify the connection between (i) 
real rotations of N-level quantum states in TV 2 — 1 di- 
mensions and (ii) their mapping to classical oscillators in 
one-dimensional physical space. Quantum evolution in 
(complex) Hilbert space is also readily compared to the 
evolution of the corresponding classical systems derived 
from Liouville and Hilbert approaches. In what follows, 
vectors are written conveniently as rows in the text, but 
are to be understood as column vectors when used in 
matrix equations. 



A. Quantum solution 

In terms of real Ai, A2 and complex V = lj\ — 1LJ2, the 
Hamiltonian for a general 2-level system can be written 
in terms of the Pauli matrices cr,; (i = 1,2,3) and cr = 1 
as 

H = ( V* Ao ) =I] W « CT «> ( 27 ) 

with wo = l/2(Ai + A 2 ) andw 3 = l/2(Ai-A 2 ). The cr 
term commutes with the other terms, so the propagator 
U(t) = e~ lHt giving the Schrodinger equation solution 
as in Eq. (|17p is readily obtained in terms of oj, cr, (i = 
1,2,3). The standard expansion of e - lU3 lTt using unit 
vector 6: = uj/lo gives 



Using the a a as the basis and inner product 
{<y a \ap) = 1 / 2Tr (cr Q &p ) = 8 af} gives ft 0a = = fl a o 
according to Eq. ©. The remaining 3x3 matrix 
giving the non-zero couplings is easily determined us- 
ing the commutation relations [crj,erj] = 2ieijk^k writ- 
ten in terms of the usual Levi-Civita tensor Cijk (equal 
to ±1 for cyclic/anticyclic permutations of the indices 
j, k, I = 1, 2, 3 and zero otherwise) summed over repeated 
indices for slightly more concise notation. Then 

Ski = (Wi^H) 

= -2e ljk (<j k \u)-<r) 
= —2u}i€ijk{crk I fz ) 

= —2l0kCijk 

(0 -U3 lo 2 \ 
lj 3 -wi . (30) 
-L0 2 Ul\ J 

The resulting equation of motion 

f = !lr = 2uxr (31) 

represents a rotation of r about uj at angular frequency 
2cj, as expected from the FVH Theorem for arbitrary 
two-level systems. The equivalence of quantum dynam- 
ics, in the case of 2-level systems, to a rotation in real, 
physical space cannot be generalized to arbitrary N-level 
systems. Representing quantum dynamics by a system 
of coupled harmonic oscillators in one physical dimension 
is general for any value of N. 
The coupling matrix is 

(-(wf + cj 2 ) WlW 2 ^1^3 \ 

wiw 2 -{ujf + CJ 2 ) CJ 2 W 3 
w 2 w 3 -(wf+wf) J 

(32) 
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giving three mutually coupled oscillators as in Fig.[TJ The real and imaginary parts of H as in Eq. (fl"8|). giving 
couplings obtained from Eq. (|16l) are 

fcij/4 = u)iL)j i ^ j 

ku /4 = L)j + LJ k - 0Ji ujj — uiiUJk i ^ j ^= k 

= ujiuj - Ui) +u} k {uj k -ui). (33) / o -u) 2 Ai wi 

For any possible ordering of the relative magnitudes of ^ = I _a } 2 I (34) 



nonzero Wj, at least one of the ku has to be negative. ^ —uil — A 2 W2 2 



C. Classical Representation (Schrodinger equation) 

The matrix f2 leading to a solution for (q, p) as a rota- 
tion e" ot of the initial state (qo,Po) is comprised of the and coupling matrix 




Four coupled oscillators are needed to represent 
(qi,q2,Pi,Pz) = {ri,r 2 ,r 3 ,r 4 ). The mutual couplings hj 
(i ^ j) given by Eq. (|16[) are the corresponding elements 
of Vi 2 . The self- couplings for i = 1,2 are 

ku = A 2 + u 2 + wf - (wi + wa)(Ai + A 2 ), (36) 

with = fcii and /C44 = k 22 . Negative couplings are 
required in general, except for the special case u 2 = 
and lvi < 0. 

The operator f2 generates simultaneous rotations in 
the planes r^-r^ associated with the nonzero flij. The 
nonzero mutual couplings represent the noncommuting 
rotations in fi. One easily shows (see Appendix) that 
noncommuting rotations share a common coordinate axis 
in their respective rotation planes, such as r 2 -r\ and 
T\-T§. Then ^21^13 = (^23) gives a nonzero mutual 
coupling k 23 . A rotation in the r±-r 2 plane does com- 
mute with a rotation in the r 3 -r 4 plane, so one expects 
the mapping from rotations to oscillators to generate at 
least one mutual coupling equal to zero. For the par- 
ticular example here, the structure of Q is such that 
(fi? 3 ) = = (^24), gi vm g zer0 f° r &13 an d k 2 4. Mass 
1 is not coupled to mass 3, and mass 2 is not coupled to 
mass 4. 

Thus, the rotation of a single vector or rotor in real 
4-dimensional space can be viewed equivalently as 4 os- 
cillators evolving in one dimension. Similarly, for N-level 
systems, the evolution of the associated rotor in N 2 — 1 
real dimensions is equivalent to the evolution of N 2 — 1 
oscillators in one physical dimension. 



D. Quantum dimer 

The quantum dimer example provided in [29j j is given 
by real V — oj\, uj 2 — 0, and Ai = A2 = ljq, giving 
uj 3 = 0. 



1. Liouville approach 

The only nonzero elements of in Eq. (f3T))) are then 
O32 = 2V = —^23, leading to diagonal entries (fl 2 ) 22 — 
(Sl 2 )33 = —AV 2 as the only nonzero elements of fi 2 in 
Eq. (|33p . Thus, two uncoupled oscillators, each with nat- 
ural frequency 2V, represent this particular quantum sys- 
tem, with the initial conditions determining the specific 
details of the time evolotion. 

For #(t) = [cx(t),c 2 (t)] and *(0) = (1,0), one easily 
obtains r(0) = (0,0,1/2) using r, = l/2 r Yt(a i p), result- 
ing in r(0) = [0, -V, 0)] from Eq. O- Then 



r(t) = - -sin2Fi , (37) 
z \ cos2Vt ) 

which is the expected rotation about axis Cj = u>i at 
angular frequency 2u>it — 2Vt given by Eq. (f^jl . Since 
the two oscillators are out of phase by 90°, the system 
can actually be represented by a single oscillator — the 
position of one oscillator automatically gives the position 
of the other from a simple phasor diagram. 
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2. Schrodinger approach 



2. Schrodinger approach 



Referring to the 2x2 block structure of f2 2 in Eq. (1331) , 
one finds off-diagonal blocks equal to zero, since they 
depend on the imaginary part of V. The two remain- 
ing nonzero blocks on the diagonal generate independent 
evolution of q and p. The q-block gives two coupled 
oscillators with mutual coupling k\2 = —2us$V and self- 
couplings fejj = (lli — V) 2 from Eq. ([33)1 . The p-block 
gives identical couplings. One can instead switch to a 
positive value for the mutual coupling, as in [29], since 
the normal-mode eigenvalues — (luq ± V) 2 of il 2 are only 
interchanged by changing the sign of V. This changes 
the sense of rotation generated by H in Hilbert space 
and hence, by fi in the real 4-dimensional space. Us- 
ing a positive coupling in this way captures the essential 
elements of the problem, but does not, strictly speak- 
ing, faithfully map the quantum system to the oscillator 
system. More importantly, negative couplings cannot be 
avoided, most generally. 

With the definitions in ^11 B 11 the initial condition 
c(0) = (1,0) corresponds to (qo,Po) = (1,0,0,0), which 
extracts the first column of f2 in the matrix multiplica- 
tion of Eq. flU to give (q ,p ) = (0, 0, —ujq, —V). The 
four oscillators must be set in motion with these initial 
positions and velocities for their displacments in a me- 
chanical implementation to correspond to the evolution 
of \*(t)) = [ci(t),c 2 (t)}. 

However, the propagator U (t) is readily obtained from 
Eq. (fT4|) in terms of the eigenvectors (1, 1) and (1, —1) for 
each 2x2 block on the diagonal, padded with zeros to 
give the appropriate four-element vector. A solution for 
the motion requires only the initial displacements. The 
given initial condition picks out the first column of U(t) 
to reproduce the solution given in [2{| . The Schrodinger 
equation requires four coupled oscillators for this partic- 
ular example, in contrast to two uncoupled oscillators for 
the Liouville representation (equivalent to a single oscil- 
lator, since they are always 90° out of phase). 



E. Symmetric unperturbed levels 



The matrix fl 2 of Eq. (|3"5j) is now diagonal for any 
general complex perturbation V. Four uncoupled oscil- 
lators, each with natural frequency (y"V w|) 1//2 , represent 
the system. Specifying | vf'(O) ) determines the initial con- 
ditions as discussed previously. This is a very simple sys- 
tem, with each oscillator evolving independently. 

The Liouville approach, by contrast, results in a rela- 
tively more complex system, albeit with one less oscilla- 
tor. Yet, for the dimer example, the Liouville implemen- 
tation is much simpler than the Schrodinger implementa- 
tion. Which approach gives the simpler set of oscillators 
and couplings depends on the specific problem. 



F. Bloch equation with relaxation 

The solution of the Bloch equation for the time depen- 
dence of nuclear magnetization in a magnetic field is rel- 
atively simple for a field along the z-axis [32[ . As is well- 
known, the transverse magnetization precesses about the 
field at the Larmor frequency while decaying exponen- 
tially at a transverse relaxation rate I/T2. The longi- 
tudinal magnetization relaxes to the equilibrium mage- 
tization at a rate 1/T\. The mapping of this motion to 
a system of damped oscillators illustrates the procedure 
described in ijllll as well as the role of non-reciprocal cou- 
plings in the model. Physical insight that may have been 
overlooked in the past is readily apparent from classical 
textbook treatments of damped oscillations. 

The inhomogeneous term F in Eq. is (0, 0, Mq/T x , 
where Mq is the equilibrium magnetization. Vector r 
represents the nuclear magnetization. Denoting ujg, as 
the Larmor frequency, the matrix L = O + R is 



r = 



1 

T 2 


-w 3 







1 

T 2 











1 

Ti 



(38) 



Consider Ai = — A2 = W3, which arises in representing 
two unequal energy levels relative to the mean energy of 
the levels. 



As described earlier, appending a column r_F to the right 
of r 2 followed by a row of zeros at the bottom gives 
Eq. (j2U)) for the oscillator equation, with 



1. Liouville approach 

There are no nonzero elements of f2 2 derived from 
Eq. (|32l) . The system is fully coupled, as illustrated in 
Fig- HI and represents the most general result for this ap- 
proach. The Schrodinger approach, discussed next, pro- 
vides a simpler representation in this case. 



L = 







V 







2lJ3 

T 2 



1 

5? 






1 





\ 





M 

If 





(39) 



and r 4 = 1 augmenting r to represent a static component 
that incorporates the inhomogeneous term TF. The nec- 
essary couplings are easily read from symmetric T$ and 
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antisymmetric Ta that sum to give T: 



( W 



Ts = 



V 









" T 2 






2uj 3 
T 2 









2 








Mo 
27 : 



2 I ; 






M 








Mo 





(40) 



A symmetric coupling k^ — — Mo/(2T 1 2 ) connected 
in parallel with antisymmetric (nonreciprocal) coupling 
734 = — Afo/(2T 1 2 ) provides the contribution to the final 
steady state magnetization through coupling to r±. 
The vanishing of k^ + 743 ensures there is no coupling 
from r3 to change the static component r^. Although 
there is no friction term in Eq. (|26[) . the mechanism that 
damps r3 is fairly transparent. Since k^i = = k^, 
r3 is only coupled to static which effectively shifts 
r 3 to z = f 3 — M , giving the equivalent equation z = 
z/T*. The self-coupling fc 33 is the source of the imaginary 
natural frequency i/T\ . resulting in the standard damped 
solution z(t) = z(0)e~ t/Tl . 

The mechanism for transverse relaxation is perhaps 
more interesting, given that the diagonal elements Ta 
(i = 1,2) cannot be the source of the damping for the 
case W3 = I/T2. More generally, since Tn — T22 for any 
value of W3, the eigenvalues are r n plus the eigenvalues 
for the antisymmetric block, which are i2iuis/T2. The 
normal mode frequencies, given by the square root of the 
eigenvalues, are U3 ± i/T2- The asymmetric coupling is 
the sole source of the imaginary frequency producing the 
required e~ t l T ' 2 decay of the transverse magnetization. 
The displacements of f\ and fi exert opposing forces on 
each other compared to the usual symmetric couplings, 
which produce oscillations, and the motion is damped. 



methodology provides the possibility for visual, mechan- 
ical insight into abstract quantum systems, as well as a 
metric for characterizing the interface between quantum 
and classical mechanics. There is a one-to-one correspon- 
dence between oscillator positions at time t and quantum 
states of the system. The formalism presented includes 
both open (dissipative) and closed systems. For systems 
represented by a density matrix, the known evolution of 
states as rotations of a single coherence vector in a real 
(but unphysical) hyperspace of A 2 — 1 dimensions [|[ has 
been mapped here to the evolution of A 2 — 1 oscillators 
in one physical dimension. The evolution of Schrodinger 
states has also been generalized here to real rotations in 
2A dimensions, which can be mapped to 2A oscillators 
in one-dimensional physical space. The scaling of quan- 
tum systems to classical systems is therefore linear for 
Schrodinger states rather than quadratic for density ma- 
trix representations. 

The results are applicable to time-independent Hamil- 
tonians, which is sufficiently general for a great many 
cases of practical interest. Time-dependent Hamiltoni- 
ans (ie, driven systems) can be approximated to a chosen 
level of accuracy by a sequence of constant Hamiltoni- 
ans over sufficiently short time steps. Modeling a time- 
dependent quantum Hamiltonian as classical requires a 
change in the spring constants of the mechanical system 
at each time step, together with a reinitialization of the 
velocities derived from the new positions of the oscillators 
according to Eq. ©. By contrast, velocities in the nat- 
ural dynamics of a system of coupled oscillators would 
not change instantaneously and discontinuously with a 
change in spring constants. 
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